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Abstract 

The maximal acceleration corrections to the Lamb shift of muonic 
hydrogen are calculated by using the relativistic Dirac wave functions. 
The correction for the 2 S — 2 P transition is ~ 0.38 meV and is higher 
than the accuracy of present QED calculations and of the expected 
accuracy of experiments in preparation. 
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1 Introduction 


This paper presents the calculation of maximal acceleration (MA) corrections 
to the Lamb shift of muonic atoms p + p~, according to the model of Caianiello 
and collaborators [|I|], 0. The view frequently held ||, || that the proper 
acceleration of a particle has an upper limit finds in this model a geometrical 
interpretation expressed by the line element 


ds 2 = g fll/ dx fJj dx u 



a 2 (x)ds 2 , 


( 1 . 1 ) 


experienced by the accelerating particle along its worldline. In (|1 . 1|) A m = 
2mc‘/T l is the proper MA of the particle of mass m, x M = d 2 x fl /ds 2 its accel¬ 
eration and ds 2 = g^dx^dx" is the metric due to a background gravitational 
field. In the absence of gravity, g^ u is replaced by the Minkowski metric 
tensor Results similar to ( |1 . 1| ) have also been obtained in the context of 
Weyl space 0. 

Eq. (0) has several implications for relativistic kinematics || , the energy 
spectrum of a uniformly accelerated particle [0, the periodic structure as a 
function of momentum in neutrino oscillations [0 , the Schwarzschild horizon 
[§], the expansion of the very early universe 0, the classical electrodynamics 
of a particle |H|] and the mass of the Higgs boson [!llj. It also makes the 


metric observer-dependent, as conjectured by Gibbons and Hawking [[12 
and leads in a natural way to hadron confinement Jl3 . 

The extreme large value that A rn takes for all known particles makes a 
direct test of Eq. (O) very difficult. Nonetheless a realistic test has also 
been suggested El- 

Using the same model, we have recently calculated JT5| in a non-relativistic 
approximation, the MA corrections to the Lamb shift of hydrogenic atoms 
and found them compatible with experimental results. In particular, the 
agreement between MA corrections and experiment is very good for the 
2 S — 2 P Lamb shift in hydrogen (~ 7 kHz) and comparable with the agree¬ 
ment of experiments with standard QED with and without two-loop correc¬ 
tions. The agreement also is good for the (l/A)Lis — (5/4)L 2 s + L 4 $ Lamb 
shift in H and comparable, in some instances, with that between experiment 
and QED (~ 30 kHz). The agreement remains good, in this instance, for D 
too. For the Lis case in D, the MA theory is worse (~ —270 kHz) than the 
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standard one in reproducing the experimental data when the two-loop cor¬ 
rections are included, but better than QED alone when these are excluded. 
Finally, the MA corrections improve the agreement between theory and ex¬ 
periment by ~ 50% for the 2 S — 2 P shift in He + . 

In this work we extend the calculation of the MA corrections to muonic 
hydrogen atoms for essentially two reasons. First, the levels of muonic hy¬ 
drogen are very sensitive to QED, recoil and proton-size effects and may lead 
to a more precise determination of the proton radius. An accurate measure¬ 
ment of the proton radius would affect all QED tests based on the hydrogen 
atom and corresponding comparisons with the MA corrections. Second, MA 
effects are larger in muonic hydrogen because the muon in the ground state is 
much closer to the proton, hence its acceleration is higher. Unlike Ref. [[EJ, 
the present calculations are fully relativistic. Section 2 contains the Dirac 
Hamiltonian, its eigenfunctions and the MA perturbations. The Lamb shifts 
are calculated in Section 3 and the conclusions are given in Section 4. 


2 The Dirac Hamiltonian 


The MA corrections due to the metric ( |1 . 1|) appear directly in the Dirac equa¬ 
tion for the muon that must now be written in covariant form (T 6 |] and referred 
to a local Minkowski frame by means of the vierbein field e M a (x). As in Ref. 


I5p and |16| one finds e a = a(x)5 a , where Latin indices refer to the locally 


inertial frame and Greek indices to a generic non-inertial frame. The co¬ 
variant matrices 7 ^(x) satisfy the anticommutation relations { 7 ^( 2 ;), 7 " (a;)} 
= 2 g^ w (x), while the covariant derivative X% = + ucontains the total 

connection u; M = \o ab u IMb , where o ab = \ [yGy 6 ], uj^ a b = (rj w e x a - d [1 e v a )e v b 
and represent the usual Christoffel symbols. For conformally flat met¬ 
rics atakes the form u ,= (l/cr)a ab 7] afl a^ b . By using the transformations 
7 M (x) = e%(x) 7 a so that 7 M (x) = <r _ 1 (x) 7 % where 7 ^ are the usual constant 
Dirac matrices, the Dirac equation can be written in the form 


ih'f (d/j, + + *Y7 M ( lnor ),M - mca(x) 


7p(x) = 0 , 


( 2 . 1 ) 


^Frorn (|2.1|) one obtains the Hamiltonian 

.3 Tic 


H = —ihea ■ V + e 7 ° 7 ^A^x) — 7 ° 7 M (lno-) )Al + mc 2 a(x) 7 0 , (2.2) 
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which is in general non-Hermitian Jl 6 j]. If cr varies slowly in time, or is time- 


independent, as in the present case, the term (Incr ) )0 can be neglected and 
Hermiticity is recovered. 

The Lamb shift corrections are calculated by means of relativistic wave 
functions [|17[ . For the electric held E(r ) = kZe/r 2 (k = l/47reo), the con¬ 


formal factor becomes c(r) = (1 — ) 1//2 , where r 0 = (kZe 2 /mAm) 1 ^ 2 ~ 

1.59 • 10~ 15 m and r > r 0 . The calculation of x p is performed classically. Ne¬ 
glecting contributions of the order 0 (Al“ 4 ) one gets a(r) r'U 1 - ( 1/2 )(r 0 /r) 4 . 
This expansion requires that in the following only those values of r be cho¬ 
sen that are above a cut-off A, such that for r > A > ro the validity of the 
expansion is preserved. The actual value of A is determined by the maxi¬ 
mum probability distance of the muon from the proton. Thus A ~ ao, where 
do = h/m^ca is the Bohr radius of the muon. The length ro has no fun¬ 
damental significance in QED and depends in general on the details of the 
acceleration mechanism. It is only the distance at which the muon would at¬ 
tain, classically, the acceleration A m irrespective of the probability of getting 
there. 

By using the expansion for a(r) in ( [2.2| ) one finds that all MA effects are 
contained in the perturbative terms 


Hr 0 = 


me 


P + ji— rja - V-7 = H + M'. 


(2,3) 


The corrections to the energy levels 2 S and 2 P are calculated by using 
the eigenfunctions of the Dirac Hamiltonian 


|^ (0) >= 


9n r k(r)Xk 

ifnA r )X-k 


(2.4) 


where are the spin functions and g nr k( r ) and f n A r ) are the radial wave 
functions 


- , ry \ 

9n r k(r ) = B k e~ p ^ 2 p 1 ~ 1 (k - iFi(--n r ; 2y + 1; p)+ (2.5) 


fnAr) = C k e- p ' 2 pi- x 


~\~n r i F\ ( n r + 1 ; 2 y + 1 ; p)] , 
Za 


iFi(-n r ;2x+ l;p)~ (2.6) 
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—n r \F\[—n r + 1; 27 + 1; p)\ . 

1 Fi(a, b , x) are the confluent hypergeometric functions. The constants in Eqs. 
( |2.5| ) and ( PTBp are defined by 


, x (A - 1 + VF ) 1 / 2 _ .. 1 AK 1 -wy /‘ 2 

B k = A(n r , k)- - , C k = -A(n r , fc)- 


A(n r , fc) = 


k — Za/XX, 

2 1 / 2 A 3 / 2 
r(2y + 1) 


k — Za/XX c ’ 




+ ?i r + 1) 


11/2 


/AA C )K)! 


The quantum numbers n r , fc and 7 are given by 

n r — n — \k\, 7 = \jk 2 — (Za) 2 , a ~ 1/137, 


(2.7) 

( 2 . 8 ) 


(2.9) 


where k is related the angular quantum number l (for instance, k — — 1 for 
the states S' and k = 1 for the states P). W is defined in terms of the energy 

E n i.j 


W = 


E n ij _ me 2 


he 


he 


1 + 


(Za) 


[n - (3 + 1 / 2 ) + (, k 2 - (Z«) 2 ]' 


- 1/2 


where n = 1, 2, 3,... , j = 1/2, 3/2,... < n, 0 < l < n — 1. Finally, 

A c = —, A = (A 7 2 - IF 2 ) 1 / 2 , p = 2Xr . 
me 

The perturbation due to Xi' vanishes, while for 7i one finds 

mc 2 r / k r°° 1 


A E = 


( 2 . 10 ) 


( 2 , 11 ) 


r 00 1 

J A -^[ 9 n r k(r ) 2 - f nrk (r) 2 }dr . ( 2 . 12 ) 


3 p + fi Lamb Shifts of the States 2Si/2,2Pi/2 


The contribution to the Lamb shift 2S — 2P is calculated by using Eq. (|2.12|) . 

For 2 S 1/2 states, one has (Z=l) n — 2, n r — 1, k — —1, and from Eqs. 
( |2.5|) and ( p. 6 |) one gets 


9 i,-i(r) = B~\e 


■p/2p7“l 



P 

27 + I 


a 

XX C 


(3.1) 
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(3.2) 


/ 1 ,_ 1 (r)=C'_ 1 e-"V 1 


1 + 


a 


P 


a 


XX C J 27 + 1 AA C 


- 2 


where the identities 18 


x 


iFi(-l, 6 ;x) = 1 - -, iFi( 0 , 6 ;x) = l 
b 


(3,3) 


have been used. Inserting Eq. (|3.1| ) and ( fh2| ) into Eq. ( [2.12| ) one obtains 
the correction to 2 S \/2 


AE(2S l/2 ) = a-! 7(0) + + 0 ^/( 2 ) 


(3,4) 


The coefficients a_i, 6 _ 1 , c_i and the integral function /(g) are defined below. 
Similarly, for the state 2Pi/ 2 one has n — 2, n r — 1, k — 1, and 


9 i, i(r) = B x e p/2 /X 1 


-1 + 


a 


P 


a 


XX C J 2q T 1 AA C 


+ 2 


/i,i(r) = C 1 e- p / 2 p ^ 1 


-1 + 


a 


P 


a 


XXJ 2 7 -T1 AA C 

Therefore, the correction to the level 2 Pi / 2 is 

AE(2Pi/ 2 ) = ai/(0) + &i/(l) + c±I(2). 

The integral function /(g) is defined as 

r 00 

/(g) = —mc 2 rtX / clpe~ p p 2 ' y ~ 2 - g , 

J2AA 

while the constant coefficients are 

P 2 ! - 


0-1 = 


b -1 = 


27 + I 


1 + 


(2 7 + l) 2 

a 


1 + 


a 

AX 


AA, 


cd 


(B 2 _ 1 -C 2 _ 1 )—-2C 2 _ 1 


c -1 = p! x 


CD 

AX 




a 1 = 


P 2 - C 2 

( 2 7 + 1) 2 


-1 + 


CD 

AX 

CD 

XX r 


+ 2 


(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 
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bi = 


2 7 + l 


ci = B\ 


-1 + 


a 


a 

AX 


(Bl - C?) 


a 


XX 


2 B\ 


c\ 


a 


(3.13) 


(3.14) 


,AA C ) 1 VAA C , 

The integral 7(g) depends strongly on the cut-off A. For A ~ ao, a numerical 
evaluation of corrections ( |3.4|) and (|3. 7| ) yields 

AE(2S l/2 ) -2.06 • 10 5 MHz , (3.15) 

AE(2P 1/2 ) - 2.99 • 10 5 MHz , (3.16) 

so that the Lamb shift correction for the muonic hydrogen atom is 

AE l = AE(2S 1/2 ) - AE(2P 1/2 ) ~ 9.3 ■ 10 4 MHz ~ 0.39meV. (3.17) 

It is interesting to note that repeating the same calculation for (p + e~) hy¬ 
drogen atoms, one finds 

AAx(p + e _ ) ~ 11.37kHz, (3.18) 

in excellent agreement with the result +10.45 kHz calculated in the non- 


relativistic approximation 15 


4 Summary 


The results of interest in the present calculation are Eqs. (+4) and (3/7). 
When A ~ ao ~ 2.6 • 10 _13 cm, the muon Bohr radius, the 2 S — 2P Lamb 
shift is given by Eq. (|3 ■ 1 7| ) . The validity of the calculation is supported by 
the value ( |3.18|) obtained for the 77-atom, which agrees well with the result 
AE L (p + e ~) ~ +10.45 kHz previously calculated using a non-relativistic ap¬ 
proximation. The result ( [3.171) is of opposite sign and much smaller than the 
Lamb shift from all sources E b = 202.070(108) meV recently calculated by 
Pachuki B, but much higher than the estimated 0.01 meV precision level 
(three-loop vacuum polarization) of his calculation. In fact it ranks higher 
than all corrections reported in |19|] with the exception of vacuum polariza¬ 
tion to leading order (205.006 meV), two-loop vacuum polarization (1.508 
meV) and muon self-energy and vacuum polarization (-0.668 meV). Mea¬ 
surements at the expected level of accuracy |1| may provide direct evidence 
for the MA corrections calculated in the present work. 


6 





















Acknowledgments 

Research supported by MURST fund 40% and 60%, DPR 382/80, the 
Natural Sciences and Engineering Research Council of Canada and NATO 
Collaborative Research Grant No. 970150. G.P. gladly acknowledges the 
continued research support of Dr. K. Denford, Dean of Science and Dr. L. 
Symes Vice President Research, University of Regina. G.L. wishes to thank 
Dr. K. Denford for his kind hospitality during a stay at the University of 
Regina. 


References 

[1] E.R. Caianiello, Lett. Nnovo Cimento 32 (1981) 65; Rivista del Nnovo 
Cimento 15 No. 4 (1992). 

[2] E.R. Caianiello, A. Feoli, M. Gasperini and G. Scarpetta, Int. J. Theor. 
Phys. 29 (1990) 131. 

[3] P. Caldirola, Lett. Nnovo Cimento 32 (1981) 264; E.R. Caianiello, S. De 
Filippo, G. Marmo and G. Vilasi, Lett. Nnovo Cimento 34 (1982) 112; 
W.R. Wood, G. Papini and Y.Q. Cai, Nnovo Cimento B104 (1989) 361, 
653, 727; B. Mashhoon, Phys. Lett. A143 (1990) 176; V. de Sabbata, 
C. Sivaram, Astrophys. Space Sci. 176 (1991) 145; Spin and Torsion 
in gravitation, World Scientific, Singapore, (1994); V.P. Frolov and N. 
Sanchez, Nncl. Phys. B349 (1991) 815; A.K. Pati: Nnovo Cimento B107 
(1992) 895; M. Gasperini, Gen. Rel. Grav. 24 (1992) 219; N. Sanchez, 
in Structure: from physics to general systems , Vol. 1, p. 118, eds. M. 
Marinaro and G. Scarpetta, World Scientific, Singapore, (1993). 

[4] H.E. Brandt, Lett. Nnovo Cimento 38 (1983) 522; in Proc. of the Fifth 
Marcel Grossmann Meeting on General Relativity p. 777, eds. D.G. Blair 
and M.J. Buckingham, World Scientific, Singapore, (1989); Found. Phys. 
21 (1991) 1285; M. Toller, Nnovo Cimento B102 (1988) 261; Int. J. 
Theor. Phys. 29 (1990) 963. 

[5] G. Papini and W.R. Wood, Phys. Lett. A170 (1992) 409; W.R. Wood 
and G. Papini, Phys. Rev. D45 (1992) 3617; Found. Phys. Lett. 6 (1993) 
409; G. Papini, Mathematica Japonica 41 (1995) 81. 


7 



[6] G. Scarpetta, Lett. Nuovo Cimento B41 (1984) 51. 


[7] E.R. Caianiello, M. Gasperini and G. Scarpetta, Nuovo Cimento 105B 
(1990) 259. 

[8] M. Gasperini and G. Scarpetta, in Proc. of the Fifth Marcel Grossmann 
Meeting on General Relativity p. 771, eds. D.G. Blair and M.J. Buck¬ 
ingham, World Scientific, Singapore, (1989). 

[9] M. Gasperini, Astro. Space Sc. 138 (1987) 387; E.R. Caianiello, M. 
Gasperini and G. Scarpetta, Classical and Quantum Gravity 8 (1991) 
659; E.R. Caianiello, A. Feoli, M. Gasperini and G. Scarpetta, in Bo- 
qoliuhovskie Ctenia , p. 134, Joint Institute Nuclear Research, Dubna, 
URSS (1994). 

[10] G. Lambiase, G. Papini. G. Scarpetta, II Nuovo Cimento 112B (1997) 
913. 

[11] G. Lambiase, G. Papini, and G. Scarpetta, Maximal Acceleration Limits 
on the Mass of the Higgs Boson, submitted to II Nuovo Cimento (1998). 
Similar results, based however on the the simple quantum mechanical 
derivation of MA rather than the more elaborate model of [1J and [^], 
have been given by S. Kuwata, II Nuovo Cimento 111B (1996) 893. 

[12] G.W. Gibbons and S.W. Hawking, Phys. Rev. D15 (1977) 2738. 

[13] E.R. Caianiello, M. Gasperini, E. Predazzi and G. Scarpetta, Phys. Lett. 
A132 (1988) 83. 

[14] G. Papini, A. Feoli and G. Scarpetta, Phys. Lett. A202 (1995) 50. 

[15] G. Lambiase, G. Papini, and G. Scarpetta, Phys. Lett. A, 1998, in press 
(|hcp-ph/9804438|) . 

[16] L. Parker, Phys. Rev. D22 (1980) 1922. 

[17] C. Itzykson and J.B. Zuber, Quantum Field Theory, McGraw-Hill Inc., 
New York, (1980); 

W.J. Grandy, Jr., Relativistic Quantum Mechanics of Leptons and 
Fields, Kluwer Academic Publisher (1990). 


8 





[18] M. Abramowitz, I.A. Stegun, Handbook of Mathematical Functions, 
Dover Publications, Inc., New York. 

[19] J.K. Pachuki, Phys. Rev. A53 (1996) 2092. 


9 



